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This study explored the relationship between urbanization, energy consumption, and CO, emission in 
the MENA countries. The panel model was utilized taking the period 1980-2009 into consideration. 
Pedroni cointegration test results showed that urbanization, energy consumption and CO2 emission 
were cointegrated. The dynamic OLS results also showed that there was a long run bi-directional 
positive relationship between urbanization, energy consumption, and CO, emission. However, the 
significance of the long run relationship between urbanization, energy consumption, and CO, emission 
varied across the countries based on their level of income and development. Moreover, long and short 
run bi-directional causal relationships were found between the variables based on the Granger 
causality test results. From the results of this study it is important for the urban planners and policy 
makers in the MENA countries to slow the rapid increase in urbanization. The level of energy 
consumption and CO, emission in the MENA countries increased more than double. Thus slowing 
down the urbanization level can help reduce the level of pollution and energy consumption. In addition, 
the increased energy efficiency, implementation of energy savings projects, energy conservation, and 
energy infrastructure outsourcing reduce the level of pollution produced by urban areas. 

© 2013 Elsevier Ltd. All rights reserved. 
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1. Introduction 


The Middle East and North African (MENA) region has 
witnessed and is still witnessing high increase in urban growth. 


* Corresponding author. Tel.: +60 174587786. 
E-mail addresses: usama81za@gmail.com, 
usama_81z@yahoo.com (U. Al-mulali), hassanhgf@gmail.com (H.G. Fereidouni), 
janicelee@utm.my (J.Y.M. Lee), csnormee@usm.my (C.N.B.C. Sab). 
1 Tel.: +60 174809810. 
2 Tel.: +60 197772188. 


1364-0321/$-see front matter © 2013 Elsevier Ltd. All rights reserved. 
http://dx.doi.org/10.1016/j.rser.2013.02.041 


In the last three decades, the urban population increased more 
than 250%. Based on the World Developing Indicators, the urban 
population in the MENA region increased more than 61% from 
total population in 2009 and this percentage is expected to rise in 
the coming years. This rapid increase in urban growth might 
cause a number of challenges such as the increase in the demand 
of energy which might be the reason behind MENA’s rapid 
increase in energy consumption and CO, emission in the last 
three decades. In addition, none of the previous studies explored 
the relationship between energy consumption, CO2 emission, and 
urbanization in the MENA countries. This motivated the researcher 
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to explore the relationship between urbanization, energy consump- 
tion, and CO, emission in the MENA region countries, namely 
Algeria, Bahrain, Djibouti, Egypt, Arab Rep, Iran, Islamic Rep, Iraq, 
Jordan, Kuwait, Lebanon, Libya, Mauritania, Morocco, Oman, Qatar, 
Saudi Arabia, Sudan, Syrian Arab Republic, Tunisia, and the United 
Arab Emirates (UAE). Exploring the relationship between the energy 
consumption, CO emission and urbanization is important because 
from this study we can find out whether urbanization is a major 
source of the energy consumption and CO2 emission increase in the 
MENA countries. 

A number of studies have explored the relationship between 
urbanization, energy consumption and CO, emission. Li et al. [1], 
Liu [2,3], Zhang and Lin [4], Wang et al. [5,6] and Parshall et al. [7] 
found that urbanization was one of the major factors that affected 
China and the United States energy consumption. Similar results 
were found in the newly industrialized countries by Hossain [8], 
in the United States by Clement and Schultz [9], in 208 developed 
and developing countries by Jorgenson [10] and in Tunisia by 
Shahbaz and Lean [11]. Donglan et al. [12] also found that the 
income and energy consumption in urban China had increased 
the level of COz emission. In addition, Poumanyvong and Kaneko 
[13] found that urbanization increases CO, emission in a number 
of developed and developing countries. While urbanization 
affects energy consumption negatively in the low income coun- 
tries it increases energy consumption in high income countries. 
The same results were obtained by Poumanyvong et al. [14] 
who found that urbanization increases transportation energy 
consumption and the relationship gets stronger the higher the 
income of the country is. Moreover, Al-mulali et al. [15] found a 
bi-directional long run relationship between urbanization, energy 
consumption and CO, emission in 84% of the investigated coun- 
tries and the relationship is affected by the level of development 
in each country. O’Neill et al. [16] found that China’s economic 
growth effect driven by the increased labor supply and associated 
with faster urbanization affected both aggregate energy con- 
sumption and CO, emission. Zhu et al. [17] revealed an inverted 
U-shape relationship between urbanization and CO% emission in 
20 emerging countries. 

Xiangyang and Guiqiu [18], Liu et al. [19,20], Li et al. [21], and 
Wang [22] found that urbanization played an important role in 
increasing CO, emission in China. Karaca et al. [23] also found 
that the urbanization in major cities such as Istanbul increases the 
pollution level in Turkey. In the Yangtze River Delta which 
represents the most developed and populated industrial area in 
China, Gu et al. [24] found that the urbanization process in this 
area increased the level of pollution. The same results were found 
in a number of developing countries by Zarzoso and Maruotti [25] 
and in China by Zhu and Peng [26]. 

York et al. [27] also found that urbanization had a significant 
impact in increasing both energy consumption and CO2 emission 
in the world. The same results were found in a number of 
developing countries by Jones [28], in China by Wei et al. [29], 
in the Asian Pacific region by Zhao and Schroeder [30], in a 
number of developed and developing countries by Parikh and 
Shukla [31] and in the Soviet Republics by York [32]. In China, 
Feng et al. [33] found that urban household energy consumption 
and CO; emission were higher than those of the rural household. 
However, Chun-sheng et al. [34] found that the CO, emission 
produced by a rural household was much larger than that of the 
urban household in China. In addition, the population, especially 
in urban areas, had a strong relationship with CO, emission; 
however, this relationship varied across countries based on the 
income level. It was negative in high income levels whereas it was 
positive in other income levels [35]. In addition, urbanization 
increased the level of energy consumption by transportation in 
China [36]. 


2. Data and methodology 


The main goal of this study is to examine the relationship 
between urbanization, energy consumption, and CO2 emission. To 
achieve this goal, the panel model will be utilized taking the 
period 1980-2009 into consideration. This model has been 
chosen because of the many advantages it has [37]. The models 
that will be built are specified as follows: 


EMi: = f(URBANit, ENCit) (1) 


ENCit = f(URBAN;+, EMit) (2) 


where URBAN is the urban population as an indicator of urbani- 
zation measured in millions of persons, EM is the total carbon 
dioxide emission from the consumption of energy measured in 
million metric tons, and ENC is the total primary energy con- 
sumption measured in quadrillion Btu. t denotes time and i 
denotes the cross section (countries). 

The urban population as an indicator was taken from the 
World Development Indicators (WDI) [38] while the total primary 
energy consumption and the total carbon dioxide emission from 
the consumption were taken from the Energy Information Admin- 
istration (EIA) [39]. 


2.1. Estimation procedure 


The first step is to examine the stationarity of the variables; 
thus, the panel unit root test will be utilized. Two types of panel 
unit root tests, namely, Im, Pesaran and Shin, and the ADF Fisher- 
type tests will be used. If the variables are stationary at the first 
difference, the Pedroni cointegration test will be used to examine 
the bi-directional long run relationship between urbanization, 
energy consumption and CO2 emission. If the models are coin- 
tegrated, the Granger causality based on the vector error- 
correction model (VECM) will be used to examine both the short 
and long run bi-directional causal relationships between the 
variables. 


2.2. Panel unit root test 


The panel unit root tests have higher power than the unit root 
tests based on the individual time series which is the reason 
behind its popularity among researchers. Two types of panel unit 
root tests, namely, ADF and PP unit root tests, will be used. Both 
tests are based on the Fisher-type tests. These tests used Fisher’s 
(1932) results to derive tests that combined the p-value from 
individual unit root tests. These tests were proposed by Maddala 
and Wu [40] and by Choi [41]. If 7; is defined as the p-value form, 
any individual unit root test for cross-section will be I. Then, 
under the null hypothesis of unit root for all N cross-sections, we 
have the following asymptotic result: 


N 
-2 X log) > Xy (3) 


i=1 


Choi also shows that 


N 
+g XL & (mi) >N(O,1) (4) 
i=1 j=1 


where Ø`! is the inverse of the standard normal cumulative 
distribution function. The null hypothesis for these tests can be 
written as 


Ho : % = 0, for alli 


Z= 
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On the other hand, the alternative hypothesis can be written as 
follows: 


F &i=0, fori=1,2,...... N; 
1) 4<0 fori=N+1,N+2,...,N 


2.3. Panel cointegration test 


The popularity of the panel cointegration grew among 
researchers due to its high power. The Pedroni cointegration test 
will be used in this study to examine the bi-directional long run 
relationship between urbanization, energy consumption and CO2 
emission. Pedroni carried out several tests for panel cointegration 
that allow for heterogeneous intercepts and trend coefficients 
across-section. The Pedroni panel cointegration follows the fol- 
lowing regression: 


Vir = GAO C+ B Xie t P2iX2it +- - -+ PmiXmi,t + ei, t (5) 


where t is the number of observations over time, and M is the 
number of independent variables. It is assumed here that the 
slope coefficients (;,..., Pmi, and the member specific intercept «; 
can vary across each cross-section. Under the null hypothesis of 
no cointegration, the residual in the equation above e;, should be 
I(1). To find out and test whether the residuals are I(1), it is 
important to run the following regression for each cross-section: 


Cit = Pilit-1 + Hit (6) 
or 
Pi 
Cir = pieit1+ X WyAeiejtVie Y) 
j=l 


Pedroni made two types of alternative hypotheses, namely, the 
homogeneous alternative (Pedroni terms it the within-dimension 
test or panel statistics test) and the heterogeneous alternative 
(between-dimension or group statistics test). The Pedroni panel 
cointegration statistics are constructed from the residuals. These 
cointegration statistics are specified as follows: 


a. Panel-p statistic 


T 5 FEON OR. 3 — & 
T\/N Zn,T-1 =rvn( yb) >, XO iini 4êi 0) (8) 


i=1t=1 


b. Panel-t statistic 

AA -1/2 N ta 
Zins = (SNT DODES DODE (9) 
c. Group-p statistic 


N/T SE T 
TN"?Ž an, r-1 =TN P Y ( és) So Cit-1 48-91) (10) 
1 f=1 


i=1 \t= 


d. Group-t statistic 


N T “We r 
“yas = a2 5 > 3 
NZ =N sp | Pa) Tae ramen) 
t=1 


i=1 


and 


The optimal mean and variance term can be applied when the 
panel cointegration test statistics are calculated. This helps the 
test statistics to be asymptotically standard normally distributed. 

If the models are cointegrated, the panel dynamic OLS (DOLS) 
will be utilized to examine whether there is a negative or a 
positive long run relationship between the variables. The panel 
DOLS was proposed by Pedroni [42]. The DOLS method involves 
expanding the cointegration regression with lags and leads 
resulting in a cointegration equation error term which is ortho- 
gonal to the entire history of the stochastic regressor innovations: 


Vix =Xit B+DIie + YS AXieyjO+Vrit (12) 
j=-4 
Adding q lags and r leads of the different regressors eliminates 
all the long run correlations between the residuals. 


2.4. Panel Granger causality 


The Granger causality test will be utilized in this study to 
examine the causal relationship between urbanization, energy 
consumption, and CO2 emission. If the models are cointegrated, 
the Granger causality based on the vector error-correction model 
(VECM) will be used. The VECM Granger causality can show both 
the short run causal relationship based on the F-statistics and the 
long run causal relationship based on the error correction term 
ect (-1). The VECM models can be presented as follows: 


l 
AURBANir = Xit + BiteCti-1 + X Sie AURBAN i1 
i=1 


l I 
+ JO Pu ACENO)t_1 + X Sit AEM) + hir (13) 


i=1 i=1 


L 
AENCit = dit + PiteCtit—1 + > it AENCie_1 
i=l 


l l 
+ $2 Pu ACURBAN) t1 + X. dix AEM) + Hi (14) 


i=1 i=1 


l 
AEM; = Hie + PiteCtit-1 + 5 Cir AEMit—1 
i=1 


I l 
+ XO Pu AURBAN) 1 + X. Su AENC)+ bi (15) 
i= i= 

A is the first difference operator, &it¢ is the constant term, Pit, it 
Qin and oj, are the parameters, ecti—ı is the lagged error 
correction term obtained from the cointegrating equation and 
Lit is the white noise error. 

On the other hand, if the models are not cointegrated, the 
Granger causality based on the Vector autoregressive (VAR) 
model which captures only the short run causal relationship will 
be used. Granger causality based on the VAR model can be 
presented as follows: 


l l 
AURBANit = dit + X Ču AURBAN 1+ X Pr AENO i1 
i=1 i=1 


l 
+ $2 ôu AEM)+ Hi (16) 


i=1 
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t l 
AENCi = it + X Gir AENCit-1+ D> Pi (URBAN) t1 


i=1 i=1 


l 
+ 5 Òir A(EM) + Hit a7 


i=1 


l l 
AEM; = it + 5 čir AEMit_1 + 5 Pir ACURBAN) it; 


i=1 i=l 


i 
+ 5 Öit A(ENC) + Hit (18) 


i=1 


3. Empirical results 


The Pedroni cointegration test requires the variables to be 
stationary at the first difference; thus the panel unit root test is 
used. Table 1 reviews the ADF-Fisher Chi-square and the PP- 
Fisher Chi-square test results which clearly show that all the 
variables are not stationary at levels. However, all the variables 
are stationary at the first difference rejecting the null hypothesis 
at 1% level of significance. This indicates that the variables contain 
a panel unit root. 

Since the variables are stationary at the first difference, the 
Pedroni cointegration test will be utilized. Table 2 reviews the 
Pedroni cointegration test results. The results gained from the 
ENC model show that ten statistics were significant rejecting the 
null hypothesis of no cointegration indicating that both urbaniza- 
tion and CO, emission have a long run relationship with energy 
consumption. Lastly, the results gained from the EM model show 
that eight statistics were significant rejecting the null hypothesis 
of no cointegration. This indicates that energy consumption and 
urbanization have a long run relationship with CO, emission. 
Thus, the Pedroni cointegration test results reveal that there is a 
bi-directional long run relationship between energy consumption, 
CO, emission and urbanization. 

Since the variables are cointegrated, the dynamic OLS (DOLS) 
will be utilized. Table 3 shows the dynamic OLS test results. From 
the results it can be seen that the significance of the relationship 
between urbanization, energy consumption, and CO emission 
varied across the countries based on their levels of income and 
development. The long run positive relationship between urba- 
nization, energy consumption and CO, emission is more signifi- 
cant in high and upper middle income countries compared to the 


Table 1 
Panel unit root test results. 


Variable ADF-Fisher Chi-square 


Level First difference 


Intercept Intercept and trend Intercept Intercept and trend 


URBAN 16.0588 28.8762 161.000° 142.618 ° 

ENC 7.24731 22.1065 111.695 ° 96.2227 ° 

EM 23.2044 19.8512 156.299° 114.394 ° 
PP-Fisher Chi-square 

URBAN 29.6784 19.3447 177.143 ° 152.318 ° 

ENC 7.90987 32.8718 394.819 ° 1016.78 ° 

EM 38.1417 34.0239 389.877 ° 906.160 ° 


Note: the unit root tests were carried out with individual trends and intercepts for 
each variable, and the optimal lag lengths were selected automatically using the 
Schwarz information criteria. 

Indicates significance at 1% and 5% levels. 


* Indicates significance at 1% and 5% levels. 


Table 2 
Pedroni cointegration test results. 


ENC model EM model 


Alternative hypothesis: common AR coefs. (within-dimension) 


Panel v-statistic 2.844874 ° 1.069360 
Panel rho-statistic — 1.695052 ” — 1.940536 ” 
Panel PP-statistic — 2.785714 ° — 5.422531 ° 
Panel ADF-statistic — 3.327934 ° — 5.422623 ° 
Panel v-statistic (Weighted statistic) 1.143215 — 1.732028 
Panel rho-statistic (Weighted statistic) — 3.236188 ° — 0.024241 
Panel PP-statistic (Weighted statistic) — 4.647079 ° — 2.646601 ° 
Panel ADF-statistic (Weighted statistic) — 4.702239 ° — 2.755977 ° 


Alternative hypothesis: individual AR coefs. (between-dimension) 


Group rho-statistic — 3.485617 ° — 3.610724 ° 
Group PP-statistic — 4.666581 ° — 5.966276 ° 
Group ADF-statistic — 4.505089 ° — 5.686275 ° 


Denotes significance at 1%, 5% and 10%; we use the automatic selection based on 
Schwarz to choose the optimal lag length. 


* Denotes significance at 1%, 5% and 10%; we use the automatic selection 
based on Schwarz to choose the optimal lag length. 

> Denotes significance at 1%, 5% and 10%; we use the automatic selection 
based on Schwarz to choose the optimal lag length. 


lower middle income countries. On the other hand, the relation- 
ship between the variables is negative in low income countries 
such as Mauritania. Thus the relationship between the variables is 
determined by the income and development levels. However, this 
study will focus on the panel effect since we are using the panel 
analysis. The panel DOLS test results show that both energy 
consumption and CO% emission have a positive long run relation- 
ship with urbanization. 1% increase in energy consumption and 
CO» emission will increase urbanization by 0.865833% and 
0.681991% respectively. In addition, the 1% increase in urbaniza- 
tion and CO, emission will increase energy consumption by 
0.571683% and 0.883922% respectively. Finally, the results show 
that 1% increase in urbanization will increase CO, emission by 
0.521975%. Moreover, 1% increase in energy consumption will 
increase CO emission by 0.826646" %. The DOLS results show 
that there is a positive bi-directional long run relationship 
between energy consumption, CO2 emission, and urbanization. 
From the results above it is clear that the increase in urbanization 
will increase both energy consumption and CO, emission in the 
long run. 

Since the variables are cointegrated, the Granger causality 
based on the VECM will be utilized. Table 4 review the 
Granger causality test results which reveal that there is a bi- 
directional causal relationship between urbanization, energy 
consumption, and CO, emission based on the error correction 
term ect (—1). The short causal relationship shows a 
bi-directional positive causal relationship between CO% emission 
and urbanization and between energy consumption and CO, 
emission. There is also a one way positive causal relationship 
from urbanization to energy consumption. 


4. Conclusion and discussion 


The rapid increase in the levels of urbanization, energy con- 
sumption, and CO2 emission in the MENA countries has motivated 
the researcher to explore the relationship between these vari- 
ables. To achieve this goal, the panel model was utilized taking 
the period 1980-2009 into consideration. The Pedroni cointegra- 
tion test results showed that urbanization, energy consumption, 
and CO2 emission were cointegrated. The dynamic OLS results 
showed that there was a positive bi-directional long run relation- 
ship between urbanization, energy consumption, and CO, 
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Table 3 
The dynamic OLS test results. 


URBAN as the depended variable 


ENC as the dependent variable 


EM as the dependent variable 


ENC EM URBAN EM URBAN ENC 
Algeria 0.393530 ° 0.699759 ° 1.468231 ° 0.852804 ° 0.560314 ° 0.608594 ° 
Bahrain 2.048767 ° 0.152598 ° 0.193759 ° 0.159218 ° 0.410222 ° 2.164659 ° 
Djibouti 2.901084 ° 2.566031 ” 2.394559 © 0.415538 € 0.085340 € 1.990326 € 
Egypt 0.110794 € 0.239291 € 1.844968 € 1.232062 € 0.736370 © 0.877338 © 
Iran 0.347081 » 0.156794 ° 0.675201 ° 2.302840 ° 0.032180 ” 0.635828 ° 
Iraq 0.310549 € 0.147413 € 0.434815 € 0.819709 € 1.293646 € 0.780348 € 
Jordan 2.325041 ° 1.205672 > 0.349752 ° 0.998112 ° 0.350243 ° 0.842465 ° 
Kuwait 0.377194 ° 0.386795 ° 0.255995 ” 0.842753 ° 0.678633 » 0.346596 ° 
Lebanon 0.869312 ° 1.206276 ° 0.429232 ° 1.081201 ° 0.498295 ° 0.874244 ° 
Libya 0.814200 ° 0.479934 P 0.593018 ° 0.524912 ° 0.135407 ° 0.837199 ° 
Mauritania —0.753241 © —0.321710 € —0.560105 » —0.976451 ° —0.215952 » — 0.999587 ° 
Morocco 0.893660 ° 0.347603 ° 0.215918 ” 1.218430 ° 0.611401 ° 0.417404 ° 
Oman 0.105560 ° 0.729784 ° 0.326549 ° 0.821088 ° 0.097488 ? 1.060345 ° 
Qatar 1.258889 > 0.408488 ° 0.319526 ° 1.241453 ° 0.602156 ° 0.984386 ° 
Saudi Arabia 0.938222 ° 0.210636 ” 0.704092 2 0.558777 ° 0.591067 ° 0.305996 ° 
Sudan 0.452048 € 0.270250 € 0.743752 € 1.346206 > 1.350971 € 1.244353 © 
Syria 1.423417 > 0.619065 € 0.119240 € 0.841189 € 1.207618 € 1.008984 € 
Tunisia 0.798178 ° 1.065079 ° 0.592585 ° 1.121626 ° 0.563196 ° 0.801184 ° 
UAE 1.291256 ° 1.866526 ° 0.282187 ° 1.402191 ° 0.222164 ° 0.699346 ° 
Yemen 0.411108 1.203529 € 0.050392 € 0.874777 ° 0.628742 € 1.052904 € 
Panel 0.865833 > 0.681991 ® 0.571683 ® 0.883922 ° 0.521975 > 0.826646 ° 


Note: the numbers presented in the table are the slope coefficients. 


* Denotes significance at 1%, 5%, and 10% levels. 
> Denotes significance at 1%, 5%, and 10% levels. 
€ Denotes significance at 1%, 5%, and 10% levels. 


Table 4 
Panel Granger causality test results. 


The independent variables 


Short run causal relationship Long run causal relationship 


AURBAN AENC AEM ect (—1) 
The dependent variable 
AURBAN - 0.252132 1.733813 © — 1.733813 © 
AENC 5.945357 ° - 1.871525 © — 3.317732 ° 
AEM 2.768967 ° 5.284918 ° - — 1.442078 ° 


Notes: the null hypothesis is that there is no causal relationship between variables, 
and A is the difference operator. ect (—1) represents the error correction term 
lagged one period. 

Indicates 1%, 5%, and 10% significant levels. 


* Indicates 1%, 5%, and 10% significant levels. 
€ Indicates 1%, 5%, and 10% significant levels. 


emission. In addition, the significance of the relationship between 
urbanization, energy consumption, and CO, emission increases 
the higher the country’s level of income and development is. 
Moreover, the Granger causality results reveal that there were 
long and short run causal relationships between the variables. The 
results clearly showed that urbanization had played an important 
role in increasing energy consumption and CO, emission in the 
MENA region since most of the energy consumption comes from 
fossil fuels which are responsible for more than 99% of total 
energy consumption. 

From the results of this study it is important for the urban 
planners and policy makers in the MENA countries to slow the 
rapid increase in urbanization. The level of energy consumption 
and CO, emission in the MENA countries increased more than 
double. Thus slowing down the urbanization level can help 
reduce the level of pollution. In addition, the increased energy 
efficiency, implementation of energy savings projects, energy 


conservation, and energy infrastructure outsourcing reduce the 
level of pollution produced by urban areas. 
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